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1 Introduction and summary 



It is common to model the boost-invariant expansion [1] of the strongly coupled (rela- 
tivistic) quark-gluon plasma (sQGP) produced in heavy ion collisions at RHIC in the 
framework of Mueller-Israel- Stewart theory (MIS) [2, 3]. One often further approx- 
imates sQGP dynamics in the relevant regime as that of a conformal theory. For a 
four-dimensional viscous conformal plasma undergoing boost-invariant expansion, MIS 
theory predicts the evolution of the energy density e(r) and the component of the 
viscous flow $(t) as (see [4]) 



3r 3r 



(1.1) 



where the two phenomeno logical parameters are: rj oc T 3 , the plasma shear viscosity, 
and Tn oc T~ x , the plasma relaxation time. 

To a large extent motivated and guided by the gauge theory/string theory corre- 
spondence of Maldacena [5, 6, 7], Baier et.al [8] and Bhattacharyya et.al [9] recently 
formulated a complete theory of the second order relativistic viscous hydrodynamics 
of conformal fluids. In this theory, the MIS linearized equations governing the boost- 
invariant expansion (1.1) are modified by the inclusion of a new term, quadratic in the 
component <£> of the viscous flow: 

3r r ' 

o -r 4r? T 4th _ Ai 2 ( L2 ) 

where Ai oc T 2 is a new second order hydrodynamic coefficient. An important observa- 
tion of [8, 9] was that for the M = 4 supersymmetric Yang-Mills plasma (SYM), and in 
fact for all (infinitely) strongly coupled conformal gauge theory plasmas, allowing for a 
dual string theory description, Ai 7^ 0. Moreover, for the near-equilibrium dynamics 1 
(where one can reasonably apply hydrodynamics at all) the nonlinear term in (1.2) is 
equally important to Tn terms of the MIS theory, introduced to restore causality in 
first-order hydrodynamics. As reported in [8, 9], and for the shear viscosity in [11], any 
(infinitely) strongly coupled four- dimensional conformal gauge theory plasma (with a 



1 The near-equilibrium dynamics corresponds to the late-time boost-invariant expansion (see [10]). 
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string theory dual) has 

77 1 2 -In 2 AiT 1 , N 

-U — , r n T = — , " J - = ^, 1-3 

S 47T 27T ?7 27T 

where s is the entropy density. 

The finite coupling corrections to rj/s and r n T for Af = 4 SYM plasma were com- 
puted in [12, 13, 14, 15, 16] and [18]: 

" 1 ^ ■ 120 — x -3/ 2 ^...\, _^_2-ln2 , 375^_ 3/2 

;i-4) 



where A is the gauge theory 't Hooft coupling. It was proposed in [19] 2 that corrections 
(1.4) are universal for all four dimensional conformal gauge theory plasmas (with equal 
a and c central charges) that allow for string theory duals. Similarly, within the same 
class of conformal plasmas, the finite coupling correction to Ai is universal. 

In this paper we report the (universal) finite coupling correction to A x . We find 

Given (1.4) and (1.5) we have now a complete set of universal phenomenological param- 
eters describing boost-invariant expansion of conformal gauge theory plasmas at finite 
't Hooft coupling. We hope these results will prove useful in numerical hydrodynamic 
simulations of RHIC (and LHC) nuclear collisions. 

The computations are quite technical, so we present only relevant steps and for the 
details refer the reader to previous work on the subject. We rely on the pioneering 
work of Janik and Peschanski [21] (and important further developments in [22, 23, 24]) 
which sets up a to study boost invariant expansion of a gauge theory plasma in a 
dual string theory setting. We use notations and results (often without quoting them 
here due to their technical nature) of [15]. Some further technical details appear in 
Appendices A and B, and supplemental data is available as [25]. 

We would like to comment on the issue of singularities in the Janik-Peschanski 
framework. To extract Ai we need to go to the third order in the late proper time 
expansion of the dual string theory background [24, 8]. This is the first time where 
the singularities of the dual background geometry can not be completely removed by 



2 The proof is given in [20]. 
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appropriately adjusting the hydrodynamic parameters [24, 26]. Specifically, appropri- 
ate hydrodynamic parameters will remove all the pole singularities in the curvature 
invariants in the late time expansion, but the logarithmic singularities will persist in 
quadratic (and higher order) Riemann tensor invariants. It was understood in [27] 
(also independently in [28]) that a singular proper time redefinition in the bulk will 
remove all the singularities identified in [24, 26]. Since such proper time redefinition 
has only a logarithmic singularity, it can not affect the condition for the absence of pole 
singularities which determines the hydrodynamic parameters. Thus, we are justified 
to use the original framework of Janik-Peschanski to extract Ai. 

Given the complexity of the computations, we feel that an independent check on 
the analysis is important. Thus, while computations in Sees. 2-4 are done in the 
Janik-Peschanski framework, in Sec. 5 we reanalyze boost-invariant expansion of a 
conformal gauge theory plasma in singularity- free approach of Kinoshita et al [28]. 
Both approaches lead to the same value of A x . 



2 Janik-Peschanski dual to a boost invariant plasma expan- 
sion 

The framework to study string theory duals of boost-invariant plasmas was proposed 
in [21]. Here we closely follow notations and analysis in [26] and [15]. Most details 3 
(including the description of the computational framework) are omitted due to their 
technical nature and the fact that they have already been explained in [15]. 

The string theory background holographically dual to a Bjorken flow of the M = 4 
plasma takes the form [21, 26, 15] 

ds 2 w = <]MNd£ M d£ N = 

(2.1) 



= e 



- 2a ^ z) g^{x)dx»dx v + e 6/5a{r ' z) (dS 5 ) 2 , 



where (dS 5 ) 2 is the line element for a 5-dimensional sphere with unit radius, and 
ds 2 = g liV dx il dx u = 



Z 2 



3 The details of the analysis are available from the authors upon request. 
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where dx\ = dx\ + dx\. The 5-form F 5 takes form 4 

F 5 = ^5 + *-F 5 , ^5 = -4 cu S 5 , (2.3) 

where u S 5 is the 5-sphere volume form. Moreover, the dilaton is = 0(r, z). 

Equations of motion for the metric warp factors and the dilaton are solved as a 
series expansion in the late proper time r — > oo, but exactly in the scaling variable [21] 

" = ^ ■ ( 2 - 4 ) 

Specifically, 

a(r,v) = (a (v) + 7a (u)^ + ^ + iai(v)J + ^ (a 2 (v) + -fa 2 (v) 

6(r, u) = f 6o(u) + 7 &o(u) J + ^ f 6i(u) + 7&i(«)) + ^ ( W + 762^) 



r 2 



(2.5) 



c(r,v) =(^c (v) + -fc (v^j + -i_^ci(u) +7Cx(w)^ + ^c 2 (v) + 7c 2 (v) 

+ ~^( C3 ^) +7c 3 (^)^ , 

a(r,u) =7ao(v) + ^ 7«i(«) + ^7«2(v) + ^a 3 (v) , 

1 1 1 . 

</>(t,v) =70o(^) + ^2/3 70i (u) + ^i73-702(w) + ^7<W) , 

where 

7 = ^C(3) («') 3 ^ ^C(3) A" 3 / 2 (2.6) 

is the leading string theory a'-correction to type IIB supergravity. As argued in 
[16], [17], the five-form is not corrected to this order. 

Solutions for the background warp factors {a^, b iy q}, i — 0, 1, 2, 3 were obtained in 
[24]; the leading string theory a'-corrections, up to the second order, {Sj, hi, &i, &i, 4>i} 
with i = 0, 1, 2, were discussed in [15]. In the next section, we extended the analysis of 
[15] to the leading string theory a'-corrections to the supergravity background at the 
third order, i.e., for a 3 , b 3 , c 3 , d 3 , 3 . 



We normalize the five-form flux so that the asymptotic AdS radius is one. 
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3 Equations of motion and solutions for {a^, 63, £3, Q13, ^3} 

We obtain equations of motion (including the constraints) at the third order for 03, 63, 
c 3 , a 3 , 03, extending the analysis in [15]. All the equations must be solved with the 
boundary conditions 



a 3 (v), b 3 (v), c 3 (v), a 3 (v), <j> 3 (v) 



= 0. (3.1) 



We find the following set of equations for the next-to-next-to-next-to-leading order 
in the late proper time expansion at order 0(7) 

,/ If,, 5w 4 -9 5w 4 -9 t . „ , nn . 

=c" + -b'l + — c' + — — 6', + 5(3 1) , 3.2 

3 2 3 (3 + v 4 )v 3 2(3 + w 4 )t; 3 V ; 

l.„ 5^ + 27 ,, , 9 + 5t; 4 A , 

° =C3 + 2° 3 + ^9) C3 + 2,(, 4 -3) a3 + 5(3 ' 2) ' (3 ' 3) 



=% + ^' + 6 » + __ ^ a' 3 + ^ ^ b' 3 + ^ c' 3 + 5(3,3) , (3.4) 



9 + 5v 4 hv & + 27~ 5v 8 + 27„, 

v{v 4 -3) a3 + v{v 8 - 9) 3 + v{v 8 - 9) " 

n -A" 4. l h» - ~ 3) n> 4. 3(^-5^-6) ^ 3(3t; 8 -l(h; 4 -21) A/ _ 72*; 2 . 
U_C3+ 2° 3 2(3 + t;> a3 + " W (t;8-9) 3 + " '" ^ 8 - 9) " ° 3 t; 8 - 9 3 
144v 2 



(3.5) 

° = 6 ' 3 + ^ + 2(, 4 -2, 2 ( +3K' 42 +2 , 2 + 3) a3 + S ^ ' ( 3 - 6 ) 
w/ 5w 8 + 27 w 32 „ _ 

=a 3 + ^8 _ 9) ^3 " + 5(3,6) , (3-7) 

0=^3 + |^|03 + 5(3,), (3.8) 

where the source terms {5(3, 1) • • • , 5(3,7)} are given in [25]. While the system (3.2)-(3.8) 
is overdetermined, we explicitly verified that it is consistent. 

Solving (3.2)-(3.8) is quite complicated. Fortunately, we do not need a complete 
solution. Our ultimate goal is to determine C from the nonsingularity of the ten dimen- 
sional metric curvature invariants to order 0(7) and to order 0(r~ 2 ) in the late proper 
time expansion. Thus we evaluate metric invariants first, find what field combinations 
affect the singularity as v — > 3 1 / 4 , and then solve just for those combinations of fields. 
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We assume that 

V = Vo + if}! + 0( 7 2 ) , C = Co + 7C1 + 0(7 2 ) , (3.9) 

and evaluate background curvature invariants to order 0(7) near 

x = 3 1/A ~v. (3.10) 

We use explicit solutions at lower orders, as well as equations of motion for the second 
and the third order to eliminate the derivatives (if possible) of 

{«3,&3,c 3 ; a 2 ,6 2 ,c 2 ,«2; a 3 ,b 3 ,c 3 ,a 3 } 
from the curvature invariants (see [15] for details). 

3.1 11 at order C(r 2 ) 

For the Ricci scalar we find 

7 f 20736(h; 13 (t; 4 - 3) 3 , 

r 2 [v 4 + 2v 2 + 3){v 4 - 2t; 2 + 3) (3 + v 4 f \ (3 + v 4 ) 4 ° 3 

- ^ 497664^3 , _ _ 

(v 4 + 2v 2 + 3)(v 4 - 2v 2 + 3)(v s - 9) (3 + v 4 ) 4 V 

+ 900*; 12 - 9882t; 8 + 7857t; 4 - 656l) a' 2 - 4 , ( 25v 28 + 324^ 

+ 189w 24 + 6696w 22 ^ + 1269w 20 + 69660i7^ 18 + 7641w 16 - 58320w 14 ^ + 22923w 12 
+ 626940r/ 2 w 10 + 34263w 8 + 542376r} 2 w 6 + 45927w 4 + 236196r/ 2 w 2 

+ 54675) In + + 2c 2 ) / _ 

y 3 + w 4 (w 8 - 9) 4 (3 + v 4 ) V 

+ 5805^5022^2187^- 1492992^ 



v 4 - 2v 2 + 3)(v 4 + 2v 2 + 3)(3 + v 4 ) 3 (v 8 - 9) 5 
x ( 37w 48 + 4434912w 24 + 14348907 + 1273320w 20 + 7452w 36 + 315801w 32 



+ 46008w 28 + 17684w 40 - 36137988w 4 - 300w 44 + 12538071v 16 - 28072332v 12 
+ 78180876^ + j* 7648 ^ ( 53^6 _ 533^2 



(v 4 - 2v 2 + 3)(v 4 + 2v 2 + 3)(v 8 - 9) 5 
+ 6978w 28 + 4518w 24 + 123228w 20 + 112428w 16 + 1289358v 12 + 1255338i; 8 



- (m279v 4 + 4153113 ) } . 

(3.11) 
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where • • • denote lower orders in the later time expansion. We recall explicit expressions 
for {a 2 ,b 2 ,c 2 } as well as present the decoupled equation for a 3 in Appendix A. Using 
(A.l), (A.4) we find 

^ = --- + ^|-3r} 3 1/4 (v / 3-18r} 2 ) + ^fj (l3VS - 4 In 2^3 - 30^ + 2C ) ^ 

- ^ 3 1/4 (4267)^ + 4C V3 - 3 - 24 In 2) 

- |7)o (75 - 78 In 2 + 876^^3 + l3C Vf) ^ 

+ y^o3 1/4 (V3 - 2 In 2^ + 84^ + C ) ^ - 64^0 lnx + finite J . 

(3.12) 

From (3.12) we find that the Ricci scalar of the string theory geometry does not have 
pole singularities as x — > + when 

^0 = ^74574' C7 = 2V31n2-i|, (3.13) 

which are precisely the conditions found from the nonsingularity of the Riemann ten- 
sor squared [23, 24], as well as higher curvature invariants [26], in the supergravity 
approximation to the string theory dual of the M = 4 SYM Bjorken flow. The differ- 
ence here (compare to [23, 24, 26]) is that f) and C are already fixed by requiring the 
nonsingularity of the Ricci scalar. 

While the pole singularities in the ten dimensional Ricci scalar are removed, given 
(3.13), the logarithmic singularity still persists. A similar observation was made at 
the supergravity level for the Riemann tensor invariants in [24, 26]. We expect that 
the remaining logarithmic singularity is removed by an appropriate change of variable 
[27, 28]. 
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3.2 n iM/p xW u,pX at order 0(r- 2 ) 

A bit more work is necessary to determine the nonsingularity condition of the Riemann 
tensor squared at order 0(7). Generalizing the notation of [26] 



+ r 4/3 



(3.14) 



+ G(r- 8 / 3 ) + 0( 7 2 ) . 

Explicitly we find: 

t [2]sugra = 1 / 1152v 5 (v 4 -3) 3 , 

3 (v 4 + 2v 2 + 3)(v 4 -2v 2 + 3)\ (3 + v 4 ) 3 ° 3 



13824v 5 t? 



(3 + v 4 ) 3 (v 8 - 9)(v 4 + 2v 2 + 3)(v 4 - 2v 2 + 3) 



8v^ 4 + 17v zu + 133v 16 + 210v 12 



+ 1674v 8 - 243v 4 + 729 ) a' 2 - "T^! ( ^ + 1^ 2 ° + 183v 16 + 60v 12 + 1647v 8 

/ K-9) 4 V 

+ 1458v 4 + 729^) (fe 2 + 2c 2 ) + ( ^^ 4 ^ 28 + 108?f v 26 + 57v 24 + 1944v 22 rf 

+ 465v 20 + 19764?f v 18 + 2565v 16 + 6480v 14 ?f + 7695v 12 + 177876?f v 10 + 12555v 8 

3-v 4 



+ 157464?fv 6 + 13851v 4 + 78732?f v 2 + 10935 j In g + ^ 4 



82944?; 12 fr / 

+ -m ^71 , P 5 w 434268v 20 + 794286v 16 + 514674v 12 

(v 8 - 9) 5 (v 4 - 2v 2 + 3)(v 4 + 2v 2 + 3) (3 + v 4 ) 2 \ 

+ 2322594v 4 - 334611v 8 - 531441 + 114642v 24 + 37740v 28 + 49v 40 + 2v 44 + 13Uv 



36 



+ 7251v 32 ) - -- -— - — 768 J VW — „ w „ -- ( v 36 - 25v 32 + 1026v 28 



(v 4 - 3)(v 4 - 2v 2 + 3)(v 4 + 2v 2 + 3)(v 8 - 9) 4 
+ 6054v 24 + 33300v 20 + 88452v 16 + 338094v 12 + 392202v 8 + 317115v 4 + 452709 



)}■ 

(3.15) 
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*W 1152^(^-3f _ 192 ^ ^ 12 

-3 ( U 4 + 2 U 2 + 3)( u 4_ 2u 2 + 3)( 3 + u 4)3 3 (3 + ^4)4 



+ 540, + 405 j a 3 - {y4 - ^ - 3)2(t;4 _ ^ 2 - 3)2(v4 - 3)3(t;8 _ Q) [ 8t; 

+ 17v 20 + 133, 16 + 210v 12 + 1674v 8 - 243v 4 + 729^) a' 2 

- , s 4 552 f (> + 183yl e + 60ul2 

(v 8 - 9) 4 (w 4 + 2v 2 + 3)(w 4 - 2v 2 + 3) V 

+ 1647, 8 + 1458, 4 + 729^ (b 2 + 2c 2 ) + - 3) ^ 

/ (3 + f 4 )° 

192u 5 Qi 



+ 



(3 + w 4)9( u 4 + 2u 2 + 3)2( w 4_2 u 2 + 3)2 3 

768t; 5 Q 2 , 

(3 + v 4 ) 7 (v s - 9f(v 4 + 2, 2 + 3) 3 (t; 4 - 2v 2 + 3)3 a2 



, 3072Q 3 ^ 8 v 

(3 + , 4 )> 8 - 9) 6 K + 2, 2 + 3) 2 (, 4 - 2, 2 + 3) 2 1 2 ° 2) 
6AQ 4 v e 3-v 4 
+ 3(3 + v A )\v* - 9) 6 (w 4 + 2v 2 + 3) 2 (, 4 - 2v 2 + 3) 2 H 3 + v 4 

i 



+ 



128Q 5 ^%5i 



+ 



(3 + v 4 ) 2 (v 8 - 9f(v 4 + 2v 2 + 3)3(v 4 - 2v 2 + 3) 3 

128^% 

3(v 8 - 9) 5 (v 4 - 2, 2 + 3) 2 (v 4 + 2, 2 + 3) 2 (3 + , 4 ) 2 
1492992v 12 Q 7 ^ 

{v 4 - 2v 2 + 3) 3 (t; 4 + 2, 2 + 3) 3 (3 + v 4 ) 5 {v s - 9) 8 
4608^ 10 Q 8 



(v A - 2v 2 + 3) 3 (v 4 + 2v 2 + 3)3(3 + v 4 ) 2 (v s - 9) 8 ' 

(3.16) 

where {Qi, • • • Qs} are given in Appendix B. Appendix B also contains explicit ex- 
pressions for a 2 , / 2 = c 2 + \i>2 [15], and the equation of motion for a 3 . As in the case 
of d 2 , although we can not explicitly solve for d 3 , we can argue that a^iy) can be cho- 
sen to be finite (along with its first derivative) as v — > 3 1 / 4 , while having a vanishing 
nonnormalizable mode as v — > + . 

We now have all the necessary ingredients to determine fji,Ci from the nonsingu- 
larity of {xf^ SUGRA (v ) + ^xf^ (v )j: using results of Appendices A and B, as well as 
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(3.13), we find 



T [2]SUGRA T [2\W 
X 3 7 X 3 



= ^J 7182 + 4 v / 23 3 /% - 155! - 72 v / 23 3 / 4 ?^ 1 

=31/4-^0+ 18 V " 

- ^7 (^V23 3/ % In 2 - 1272v / 23 3/4 In 2 + SOOgOv 7 ^ 3 / 4 + 2iVsS 2 ~ 21v / 23 3/ % 

\ i ql/4 / 

+ 83 1/4 >/2Ci - 176^1 + 32 v / 23 1/ % — + — - 7 245 2 + 1104^ + 16>/23 3/4 Ci 

/ x 4 72 \ 

- 51v^3 1/4 5i - 7632 v / 23 1/4 ln2 + 144630v / 23 1/4 + 64v / 23 3/4 <5 3 + 24 v / 23 1/4 <5 1 In 2^ \ 

+ J_ 7 ( _4^ 2 + 264^ + 2 v / 23 3/4 Ci + 22 v / 23 1/ % - 954v / 23 1/4 In 2 + 22896v / 23 1/4 
36 \ 

+ 8V23 3 / 4 5 3 + 3y/23 1 /% In 2^1 4 - —7 ( -9V23 3/ % In 2 - U10V23 3/4 In 2 

/ x l 36 \ 

+ 30384>/23 3/4 + 44^2 + 24 v / 23 3/4 5i + 6>/23 1/4 Ci + 264>/3t7i + 24v / 23 1/4 <5 3 ' 



+ ^8-\/23 3/4 - ^ (-3V23 3 '% - 49975v / 23 3/4 - 216^1 + 12-s/3<5 2 ) 7^ lnx + finite . 

(3.17) 
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(3.18) 



The residues of all pole singularities in (3.17) vanish provided 

= 3VV2 . _ 25/2 3/4 , 

432 v ; 

Ci = - V ^ 33/4 f 122238v / 23 3/4 + 18V23 3/4 5i In 2 - 67v / 23 3/ % + 64^2 
216 \ 

- 5724 v / 23 3/4 ln2 + 144 v / 23 1/4 5 3 ^ , 
<Si = -288, 5 2 = -14W23 1 / 4 , 

where we kept explicit dependence on Si , <5 2 in fji and Ci (as determined by the van- 
ishing of the residues of the poles in (3.17) up to order three inclusive); the vanishing 
of the residues of the second order and the first order poles in (3.17) determines Si and 
S 2 - 

As for the Ricci scalar (3.12), the logarithmic singularity in the Riemann tensor 
squared at the third order in the late-time expansion, (3.17), remains. This remaining 



5 The value of r)i agrees with the one determined from the nonsingularity of the second order 
late-time curvature invariants in [15]. 
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singularity has both the supergravity piece (in agreement with [24]) and the new 0{j) 
contribution. 

Notice that while fjo, Co (3.13) are determined unambiguously from the nonsingular- 
ity condition of the background geometry, the absence of singularities is not a powerful 
enough constraint to fix C\ 6 . This fact will not preclude us from computing a definite 
value of non-linear second-order hydrodynamic coefficient Ai. 

3.3 n^'RT at order 0[r- 2 ) 

Analysis of the square of the Ricci tensor can be performed in the same way as for the 
Riemann tensor squared. We find 

n^TV" = ■■■ + ^j^i^v^ lux + finite - 1920 d 3 (a;)| , x = 3 1/4 - v -> 0+ , 

(3.19) 

where we explicitly indicated the dependence on order three fields; as before, • • • in- 
dicate lower orders in the late time expansion studied in [15]. We pointed out above 
that 0:3(1;) can be chosen to be finite as v — > 3^ 4 ; this would guarantee the absence of 
pole singularities in 1Z I1V 1Z' 1U to orders 0(t~ 2 ) and 0(7). 

3.3.1 Higher order curvature invariants 
As in [26] we denote 

^ \ivp\ = ^ fiiuifiu ■ ^ 1X1 1 px 5 (3.20) 

where 

n [0] ^ pX = n^ pX . (3.21) 

We further define higher curvature invariants Z' 2 "', generalizing (3.14): 

= (it ]sugr \v) + 7 4 2nlw ») + ^ (ir ]suGRA (v) + 7 4 2 " iw ») 
+ ^ (ip sugra (v) + 7 4 2ni %)) + y 2 (ip sugra (v) + 7 ir ]w (v)) 

+ 0(r- 8 / 3 ) + 0( 7 2 ) . 

(3.22) 



'Although, unlike the analysis up to the second order [15], 171 is determined unambiguously here. 
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Given the complexity of the analysis, we checked at order 0(r 2 ) only the nonsin- 
gularity of X^. Using the results of the Appendices A and B, as well as (3.13) and 
(3.18), we find 



T [2\SUGRA T [2]W 
X 3 t- 7-i-3 



/640 5636512 r-\ , 

(^T7lV2 + 7 35/4 V2j In x + finite. (3.23) 



XEE3 1 / 4 — u->0 + 

The supergravity part of the logarithmic singularity in (3.23) agrees with the corre- 
sponding computation in [26]. 



4 Ai for the Bjorken flow of J\f = 4 SYM plasma 

In the previous section we analytically evaluated a'-corrected supergravity background 
dual to the Bjorken flow of Af = 4 SYM plasma at finite coupling to order 0(t~ 2 ) in 
the late proper time expansion. We can now extract the boundary energy density e(r) 
from the one-point correlation function of the boundary stress energy tensor using the 
a'-corrected holographic renormalization developed in [29]. We confirmed that the final 
expression for the energy density can be evaluated as in the supergravity approximation 
[21, 23, 26]: 

N2 v 2a(u,r) „ 1N 

£(t) = o hm j \J . (4.1) 

Using the results of the lower orders in the late proper time expansion [15], the details 
presented in Appendices A and B, (2.5), (3.13) and (3.18) we find 7 

A 2 (6 + 576 7 + 7 6 X ) 1 A^ 2 2 1 / 2 3 1 / 4 (15667 + 8 + 75i) 1 



e(r) 



12tt 2 t 4 /3 48tt 2 



T 



2 



+ gg^j- ( 12 + 24 In 2 + 7 (2^ In 2 + S ± + 7086 + 4212 In 2) J + C(r~ 10 / 3 ) . 

(4.2) 

The string theory result (4.2) should now be interpreted within second order rel- 
ativists conformal hydrodynamics [8, 9]. For the Bjorken flow of the M — 4 SYM 
plasma we expect [8] 



c gauge 



J) 



~ 4 / 3 - 2 Vo r- 2 + 



\vl - I (Vor* - A?) 



C 

3/ 4 n /e\-i/4 . ., n /e\V2 



8/3 + ^ ( T -10/3) ? 

(4.3) 



7 Dcspite the fact that at order three 5i and 62 are fixed (see (3.18)), we keep them arbitrary to 
compare with [15]. As in [15], the dependence on Si disappears in physical quantities. 
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where C is an arbitrary scale, related to the initial energy density of the expanding 
plasma. 

To match the string theory result (4.2) with (4.3) we need to recall the equation of 
state for the M = 4 SYM plasma [30] 

e (T) = K 2 N 2 T A (1 + 15 7 ) , (4.4) 

o 

and the M = 4 SYM relaxation time r n [18] 



_ 2-ln2 375 , A . 



Ultimately, we find: 



2 = J_ (i + 1207) , ^ = J- (l + 215 7 ) . (4.6) 
s An Tj Zty 

Notice that the ratio of shear viscosity to the entropy density agrees with the results 
reported in [12, 13, 14, 15], and the supergravity part of Ai agrees with computations 
in [8, 9]. 



5 Computation in the framework of Kinoshita et al 

In this section we compute the finite coupling correction to Ai by using the framework 
of Kinoshita and collaborators [28] for finding the holographic dual of an expanding 
plasma (see also [27]). We will work with the five-dimensional action 

S = ^J<fxy/=g(R + 12 + TW), (5.1) 

W = ~Ca^C^ ef C g ^C^ + C^% c f Cj' h e C^ , (5.2) 

where C a b c d is the five-dimensional Weyl tensor. The use of this action was justified 
in [20], where it was also shown that it leads to universal finite coupling corrections 
to hydrodynamic coefficients. The holographic dual to the Bjorken flow of the CFT 
plasma is taken to be [28] of the form 

= - r \dT 2 + 2drdr + r Ve" (l + - V dy 2 + r 2 e c dx\ . (5.3) 
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Plugging this form for the metric into the Einstein equations of motion, one can find 
the functions a(t,r),b(t,r),c(t,r) order by order in a late-time expansion as before: 

a(t,r) = a (v ) + ua\{v ) + u 2 a 2 (v ) + ... , 
b(t,r) = b (v) + ub 1 (v) + u 2 b 2 (v) + ... , 
c(t,r) = c (v) + uci(v) + u 2 c 2 (v) + ... , 

where v = rt 1//3 ,w = t~ 2 ^ 3 . In [28] these functions were explicitly computed: 

4 

W 

a (v) = l- — , bo(v) = cq(v) = 0; (5.4) 



v 



2 ^ 1 + l)v 4 -3M 1 w 4 v + w% 6 + 1 

W ) = -q 5 ' M" = , 

26 1 / w 4 \ arctan (^) + ±m(^) -f l0 ' 0J 

-ID ^ >,ln 1 — • ' 



3v 2 V v 4 7 3w 

2(t; 4 + W 4 )e 2 4 (u 3 -3iu 4 Mi) 6 £ 2 (^ 4 - 3u> 4 ) 
^ " 3^ 9^ + 9^ 

| (v 4 + w A ) ( 9w 2 M 2 + 1) arctan (£) ( w 4 - 2u> 3 v + u> 4 ) (9ty 2 M^ - 1) ln(w - w) 



6t> 5 u> 12v 5 w 
(w 4 + 2w 3 v + w 4 ) (9w 2 M 2 - 1) \n(v + w) (9M 2 w 4 + w 2 ) In (v 2 + w 2 ) 

+ 12v 5 w + o^ 4 

3(3(12 \n(v) + 5)vAfi + 4) AW + 4 (3wM 2 w 4 + v 3 ) 

18^ ' 

i 2 1 / , N 4 tt\ (9w 2 77 2 - 2w Mi + 1) arctan (±) 

+ 
+ 



Qv 2 4 \ v w J 4vw 

{3wM l - l)(2v -3w + 3w{4v - 3w)M 1 ) \n(v - w) 

2Avw 2 

(3wMx + l)(-2u - 3w + 3w(4:V + 3w)M l ) \n(v + w) 

2Avw 2 
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, (12wM lV 5 - 6wv 4 + ir (v 4 - w 4 ) v + 2w 5 ) M lW 3 2£ 2 

C ^ V ' ~ 3 ( v s _ vw *f + 9^ 

4t; 2 ln(i;)Afi (6 (w 4 - hv 4 ) M lW 4 + 4v 3 (v 4 + w 4 )) 2£ 2 
+ 3 {v 4 - w 4 ) + g ( w 5 _ vw 4^ +3^ 

(v 4 + 3w 4 - 3w 2 M 1 (Avw 2 + 9 (v 4 - w 4 ) Mi)) arctan (£) 
+ 18v 2 (v 4 -w 4 )w 

((3wM 1 - 1) ((v + w) (v 2 - 2wv + 3w 2 ) - 9(v - w)w (v 2 + w 2 ) Mi)) ln(w - u?) 

36f 2 (f — u>) (f 2 + w 2 ) w 
((3wMi + 1) ((u - u>) (w 2 + 2wv + 3w 2 ) + 9w(v + w) (v 2 + w 2 ) Mi)) \n(v + w) 

36v 2 (v + w) (v 2 + w 2 ) w 
(3Miv 3 + w 2 ) In (v 2 + w 2 ) ttv 3 - 3w {4M 2 w 4 + v 2 ) 
9v 5 — 9vw 4 9 (v 5 — v w 4 ) w 

(5.6) 

Both £i and £ 2 are gauge degrees of freedom which can be set to a convenient value. In 
what follows we take £i = — 1 and leave £2 arbitrary as a cross-check on our calculations. 

As explained in [28], the energy momentum tensor of the plasma can be read off 
from the function a(t,r), by expanding the function a(t,r) in the large r limit. More 
concretely, 

a{t,r) = 1 + ^ + ... + y + 

3 N 2 

Trr = ^) = -ii-^ (4 V)- ( 5 - 7 ) 

Performing this expansion in the above solution to a(t,r) leads to 
3iV 2 



(^-2M^ + p + 4M2 )rj. (5.8) 
This is to be equated with the hydrodynamic expansion for e(r), 



< T ) _ ,.-4/3 ,„ T -2 , H 2 + 4(A°- % 7-») __ 8/3 

— r — z?] r H r 

e 6 



(5.9) 



e = -—uj 4 , Mi =770, M 2 = A°-77or°. (5.10) 



leading to the identifications 

3N 2 
2 4tt 2 

These constants are in turn fixed by imposing regularity of the function c(t, r) [28]. An 
expansion of this function around v — w has a simple pole, unless we set 

1 ln2-l . 
Mi = — , M 2 = — . 5.11 

3w 6w 2 
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Finally, the hydro dynamic coefficients are obtained by the relations 

/ \ 3/4 / X -1/4 / f \l/2 

r ] = e r ]0 (-\ , rn = r°f-J , Ai = e A? ( -J . (5.12) 

We now want to compute 7 corrections to these by using the 7 corrected action 
(5.2). Since it is not practical to find the equations of motions from the action (5.2), 
we will work in an effective action framework, as in [15]. To do this, first notice that 
in the metric (5.3) there are two implicit constraints, namely 

g rr = 0, # rr = l. (5.13) 

It is not correct to impose these constraints at the level of the action. They should 
only be imposed on the equations of motion. Therefore we modify the metric (5.3) to: 

ds 2 = g^dx^dx" 



-r 2 a(t, r)dr 2 + 2(1 + h(t, r) - -a(t, r)g(t, r) ) drdr + av ' 2 ' dr 



+r 2 r 2 e 2b(t,r)-2c(t,r) + ^ dy 2 + rV(dx? + dx\) . (5.14) 

We further substitute 

a(t,r) = ao(v) + ua\(v) + u 2 a 2 (v ) + "ya(t, r) , 
b(t,r) = b (v) + ubi(v) + u 2 b 2 {v) + jb(t,r) , 
c(t, r) = co(v) + uci(v) + u 2 c 2 (v) + 7c(t, r , 

and evaluate the action (5.2) on the modified metric. A few comments are in order: 

• One only needs to evaluate the action to linear order in g, h. 

• It is sufficient to compute W to linear order with respect to each and every field 
a, b, c, g, h, meaning no mixed terms such as ab, gc can appear in W. 

• Since we are interested in computing a, b, c to quadratic order in u, it is sufficient 
to evaluate W to linear order in a 2 , b 2 , c 2 and quadratic order in a±, b±, c\. 

Variation of the action S = S(a, b, c, g, h) with respect to the various fields lead to 
the equations of motion 

5S 5S 5S . 

0, — =0, — — - =0, (5.15) 



5a(t,r) lg!h=Q ' 5b(t,r) lgh=0 ' Sc(t, r) lgh=0 
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and to the constraints 

=0 , =0. (5.16) 

Sg(t,r) lgA=0 5h{t,r) lgA=Q 

After finding these equations, one expands a, b, c to quadratic order in u, exactly as it 
was done for the unhatted quantities. This leads to a set of equations for the hatted 
quantities hi, b\, Ci, i — 0,1, 2. 

5.1 Order - 

It will be convenient to perform the change of variable y = w/v. In terms of this 
coordinate the equations of motion are 

in oO 
°o — D B ' 



a;/ _K + 12«o_ 2(^3)c; +2(1 _^ = ^ (517) 

y y y 

The solution of this system of equations is straightforward. Imposing the boundary 
conditions a (0) = &o(0) = c (0) = leads to 

ao(y) = ^?y 4 + + ^y 12 (72 - 49y 4 ) - (l + y 4 ) , 

bo(y) = -^y 12 - B°y , 

co(y) = -^-^ 12 + C 1 °ln(l-y 4 ). (5.18) 

The coefficient B® is a gauge parameter [28] and we will set it to zero in what follows. 
The constraint equations further impose A^ = = 0. Quite generally, the constraint 
equations can only affect the coefficients of the solutions to the homogeneous equations, 
and since these come from the supergravity part of the action, C 4 should not alter the 
constraints from those at the supergravity level. Another way to see this is to note 
that one can always expand the equations of motion about y — 0, where all the source 
terms are negligible. The constraint equations must necessarily be the same around 
y = as at any other point, as they only fix constants. 

Expanding ao(y) around y = and taking the coefficient of y A we see that this 
solution modifies eo at order 7 by some undetermined constant A\, which will however 
not affect any of the physical results as we will see. 
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5.2 Order - 1 



The equations of motion are exactly the same as before (5.17), performing the replace- 
ment — > 1 in the indices of functions and sources. The solutions are now: 



W a M =Atf + AW - 10(-343^-W + 360, + 216)^ 

2(Aly 5 + B\(l + y*)y) 
3 

720y 12 (2y + l) 



wb 1 (y)=-B 1 1 y + 
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A , , 960y 13 2873y 12 490y 9 335y 8 52y 5 555y 4 

k; ci (y) = 1 1 1 

yy> 11 66 99 22 11 22 



+ ^(-528^ - 2160)y + ^ (qQ(A\ + 15) arctan - 33(A? + 15)tt) 

(llA° + 34%(y 2 + y + l) 

66 (y 3 + y 2 + y + 1) 
(15 + A\) arctanh(y) + (2A\ - 3 (80 + C{)) In (1 - y 4 ) 

12 ' 

(5.19) 

The constraint equations now imply A\ = A\ + C\ = 0. We may further set B\ to zero 
as we did at order zero. The constant C\ is determined by imposing regularity at y — 1 
(that is, v = w). The function c\(y) has a pole there unless we set C\ = (A® — 165)/2. 

Expanding a\(y) around y = and taking the coefficient of y 4 we get a non-zero 
contribution. This implies that to this order the energy density is now given by 

3 ^...4^1 ,0-,--4/3_ O ^ /f ^ ^-165 \ 2 



c(r) = - ^ 4 ^(1 - Ah)r-*« - 2 (M x + j r 2 j , (5.20) 

with Mi given by (5.11). Now, in the late time regime we assume that the expressions 
for the energy and entropy densities are given by 

e = -7r 2 iV 2 T 4 (l + 15 7 ) , s = — iV 2 T 3 (l + 15 7 ) . (5.21) 
8 2 

Then, this result together with (5.12) leads to 

1= 1 (1 + 1207), (5.22) 

S 47T 

in perfect agreement with what was previously known in the literature [12, 13, 14, 15]. 
Notice the undetermined constant A\ has canceled out of the result. 
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5.3 Order - 2 



The equations of motion are as in the previous sections, but now using the source terms 
of order 2. These are too long to display here, but can be obtained from the authors 
at request. For small y all source terms vanish and so we get only the homogeneous 
equations solutions, namely 

w 2 a 2 (y) = A& + A 2 2 y\ 
w 2 b 2 (y) = -B 2 y, 
w 3 c' 2 (y) = C 2 y 5 . (5.23) 

One may once again choose B\ = 0. The constraint equations force A\ = 2Af + Cf = 0. 
Therefore to determine A 2 , which is the parameter that enters the energy density 
expansion, we must find out the value of C\. This can be determined by imposing 
regularity of the full solution at y — 1. Close to y — 1 one has 

, (1 + 2 In 2)^ + 18(7? -510 In 2 - 1085 ,„ , 
c' 2 (y) = K - ' 1 72 ^ _ y) + (Regular at y = 1) . (5.24) 

Imposing regularity forces us to pick 

d = 1085 + 5101, 1 2-.4; ( l + 1 n4) = 
18 

The energy density expansion coefficient of r~ 8 / 3 receives a modification, namely 

9 V 2 + 4(A? - r^or") = 9M 2 + 4 M 2 
6 6 

Then, using (5.12) and the known result [18] 



- 1 A{. (5.26) 



m 2-ln2 375 . „. 

rnT=— + - 7 , (5.27) 

together with the expression for rj/s found in the previous section, leads to 

^ = ^(1 + 215 7 ), (5.28) 

in agreement with (4.6). 
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A {a 2 , &2, c 2 } and a 3 

Explicit analytic solutions for {02,62^2} were found in [24]: 

_ (9 + 5v 4 )v 2 (9 + t;> 4 „ 2 (-1053 - 171v 4 + 9v 8 + 7v 12 )v 4 
tt2 ~ 12(9 -v*) ~ 72(9 -v 8 ) +V 6(9 - v*) 2 



1 , VZ-v 2 3„ 2 3 
H — In — — -77 In 



v 



4 



4 



8^ VS + v 2 4' 3 + v 
tt 2 v 2 (9 + v 4 ) „ v A „ 2 (-9 + 54v A + 7v 8 )v 4 



288^3 12(9 -v*) 72(3 + v*) 1 6{3 + v*){9 - v*) 
1 



li 2 




(A.l) 



4(3 + v 4 ) 



2 



(39 + 7v 4 )v A 



+ f} 2 




+ ~7V l n 



4(3 + v 4 ) 24(3 + v 4 ) 1 2(3 + v 4 ) 2 8^3 V3 + v 2 
3-v 4 



4'' ~"3 + v A ' 
where {fj, C} are arbitrary parameters. 

At third order (for the warp factors {03,^3,03}) one obtains four second order 
ODE's and a first order constraint, all linear in {03, 63, C3} (and their derivatives). The 
two additional equations are the constraints used to fix the radial coordinate and the 
late time in the boost-invariant metric ansatz (2.2). It is straightforward to use these 
constraints to solve algebraically for b' 3 and c' 3 in terms of a' 3 and the lower order warp 
factors. These expressions can further be used to obtain a decoupled second order 
ODE for a 3 : 

„ 5v w + 18v 12 + 2l6v 8 + 126v 4 + 243 , 
° _ ° 3 + ( V s-9)v(v* + 2v 2 + 3)(v 4 -2v 2 + 3) a3 + Jaa ' ( } 
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2f/(5v 12 + 9v 8 + 27v 4 + 135) 3 - v 4 

aa ~3{v 4 - 3) 2 {v 4 + 2v 2 + 3){v 4 - 2v 2 + 3) H JT^ 4 

32V3v 6 f)(v 8 + 18v 4 + 9) v 2 
+ {v s - 9) 2 {v 4 + 2v 2 + 3){v 4 - 2v 2 + 3) aiC 71 

- 7tt- a - rj-r^ o 77-5 \T (295245 + 393660?/ - 53k; 32 

3{v 4 - 2v 2 + 3){v 4 + 2v 2 + 3){v s - 9) 5 \ 

- 255Qv 28 - 9882?/ 4 - 21708?/° + 27702?/ 6 + 358668?/ 2 + 898857?/ + v 40 

- 42768r} V 6 - 1113264?/ 2 ?/ 2 - 2554416?? 2 ?/ 8 - 7336656?/ 4 ?) 2 - 12912048?/ 2 / 

- 46294416?/ 2 ?/ - 9447840?? 2 ?/ + 432?/°C + 12v 34 C - 52488t/ C + 27216?/ 4 C 

- 78732?/C - 3024?/ 2 C , + 1296?/V + lUv 34 fj 2 - 314928tj 6 C + 648u 26 C^ . 

(A.3) 

For the computation of C we actually need only the asymptotic solution of a' 3 (v) as 
x = 3 1 / 4 - v -> 0+. Using (A.2) and (A.3) we find: 



da 3 (v) 



dv 



f, 3 3 3 ' 4 1 / Uf, 3 V3 rj V3f)C\ 1 1 _ ± 

—5— — + I —ir, — + -T7T- I — + ^3 — + 0(x ), 

(A.4) 



! = i , ,_ ( ,_ ()| 8 x 4 \ 24 12 144 J ./•• x 2 



where A3 must be fixed so it satisfies the boundary condition (3.1). 

B Data for computing 

The coefficients Qi in (3.16) are given by: 

Qi = 102036672 + 1771475i + (-177147<$i - 102036672)?/ + (-4527095i 

- 8039875644)?/ + (-229635^ + 17670478860)/ 2 + (4374^ - 10087292556)/ 6 

+ (2301024348 + 4228250?/ + (14094^ - 758889972)?/ 4 + (373696740 + 1625i)t/ 8 

- (9455i + 73262340)?/ 2 + (-2075i + 3556980)?/ 6 - 9?/% + ?/% , 

(B.l) 
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Q 2 = -258280326t) 5i - 74384733888^ + 430467215 2 + (-86 9 3 4 427/ 5i 
+ 413248521607}o + 526126595 2 )tj 4 + (-1004423490r/ 5i + 9228872608632r) 
+ 113196933<5 2 > 8 + (1581922715 2 - 1491223446^1 - 33853663263240r} )t; 12 
+ (1053237335 2 - 798224382^^ + 69144574698756r} )t; 16 - (388879096 14036?7o 

- 324835115 2 + 170809074r} 5i)^ 20 - (14604786r7o5i273375<5 2 

- 22636564766388r} )^ 24 - (47042375 2 + 12251574r} 5i + 5034336682644r} )v 28 

- (22295255 2 + 5466042?^! - 4161674125224r} )t; 32 - (5682155 2 + 906940950984r} 

- 1465290r/ 5i)v 36 + (371486277216t} - 765455 2 + 1625994r} 5i)v 40 - (49680556848^ 

- 68135 2 - 451278^iK 4 + (47142t} 5i + 18898470612r} + 93995 2 )v 48 + (38215 2 

- lllO^ - 2382365412r} )w 52 + (-342r} 5i + 9235 2 + 258200244r} )w 56 + (129<5 2 

- 4642164r}o + 126r} 5i)w 60 + (85 2 + 24^5^ , 

(B.2) 

Q 3 = -387420489r} 5i - 148769467776^ + 43046721<5 2 + (-1607077584r} 5i 
+ 1530550085 2 - 661197634560%)^ 4 + (2210707403676^ - 3386342052^0^ 
+ 2678462645 2 )^ 8 + (10723033860480r} + 2324168645 2 - 3409725456?^iK 2 
+ (-2092224168^1 + 1030601885 2 + 18352430617176r} )^ 16 + (186857285 2 

- 6136825576320r}o - 970923024r} 5i)w 20 + (15039362941668r} - 4846392<5 2 

- 352115748?}o5i)t; 24 + (-1235549044224^ - 51554885 2 - 74427984t} 5i)w 28 
+ (2004210714960^0 + 6647022r/ 5i - 22156745 2 )^ 32 + (11706768?7o5i 

- 131529982464^0 - 572832<5 2 )w 36 + (4706100r/ 5i + 188278471332?7o - 598325 2 )v 40 
+ (256325 2 - 7695285120?7o + 1178064?7o5i)v 44 + (157085 2 + 224640?7o5i 

+ 2953736280r7o)^ 48 + (34128^1 + 39365 2 + 208054656?7o)t; 52 + (3348^^1 + 504<5 2 
+ 6959196^ )^ 56 + (144r7o5i + 325 2 )v 60 + (5 2 + 3?7o5i)^ 64 , 

(B.3) 



23 



Q 4 = (6025163444928?$ + 69735688020?$<5i - 13947137604?$<5 2 )?j 2 - 278942752080t) 
+ (539289320688^5i - 86782 189536^5 2 - 301239576062928^)w 10 - 6457008155 2 
+ 3874204890t} 5i + (378122397264t) - 1822311189<5 2 + 12957063021t) 5iK 
+ (34044733308852t} - 23681010965 2 + 296352759247} 5i)^ 12 + (-163296r$<5 2 

- 869 1 66 2 887$ - 34992?$5i)v 58 + (518856477120?$6i - 75303063936?$5 2 

- 1212305862965760?$)?/ 4 + (6585380676t) 5i - 63131883819804t} - 83613384<5 2 )?j 24 
+ (-1405753115292t}o - 73136196t) 5i + 21636725 2 )t/° + (-7464960?$ + 5184?$<5 1 

- 10368?$5 2 )?/ 2 - (6931186200t}o + 762845 2 + 10765 447} 5i)w 52 - (49589822592?$<5 2 

- 2727440242567$^ - 35704672266240??o> 6 + (-15014979725 2 - 122924025168840t} 
+ 256941094687} 5 1 )tj 16 + (4697466048?$<5 1 + 16703781 12?$<5 2 + 134641969016832?$) 

x v 30 - (5<5 2 + 15t7o^i + 151207}o)7J 68 - (1275264?$5 2 + 23388092928?$ + 1135296?$<5i) 
x v 54 + (-5089392?$5 2 - 322667880480?$ - 14082336?$5i)t; 50 + (292743211968?$<5i 

- 33391500912?$<5 2 - 200610438731 1648?$y 8 + (-127£ 2 - 384048t/ - 5227} 5i)t; 64 
+ (-58217292t} - 10116?yVi - 16085 2 )tj 60 + (-6054175872?$5 2 + 116968038336^1 
+ 658418155614720?$)t; 22 + (446541665 2 + 183812976t/ 5i - 16291341648432?7;o)?/ 2 
+ (156045249367} 5i - 5968672925 2 - 1199864659094647}o)w 20 + (19385568?$5 2 

- 20320117294896?$ - 547029936?$<5i)?/ 2 + (185597568?$5 2 - 1635526080?$5i 
+ 14338868355072?$)t; 38 + (-6316380??o5i - 2729165 2 - 58737530376?7/o)?/ 8 

- (26068068t/ 5i + 275114490228t/ + 3440885 2 )w 44 - (110621376^1 + 8304768?$<5 2 

- 825111360000?7o> 46 + (-36868394133108r/ + 1799419860t/ 5i + 584191445 2 )w 28 
+ (148847225 2 - 1 05925 158?yVi - 5478599707344r/ )w 36 + (-2153635128?$6i 

+ 717878376?/^2 - 215937884784960?$)?; 34 + (348880387687$^ 

- 1623120631374384?$ + 1570231008?$5 2 )?j 26 + (-233280?$ - 324?$5 2 + 324t$5 1 )tj 66 
+ (107373276t}o - 133685 2 - 127188?//o5i)t; 56 + (-25636713845 2 + 23819185620t/ 5i 

- 8131232879197277o)t; 8 , 

(B.4) 
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Q 5 = 5782609521v 4 + 3314566899v 24 + 6632626437v 20 + 10567349991v 16 
+ 12018183921v 12 + 9093486951v 8 + 1375999893v 28 + 26344593252r? 2 v 2 

- 146014477632?? 2 v 6 + 10005971 148?? 2 v 10 - 59066478504v 14 r? 2 - 75754434492?? 2 v 18 

- 55948533840v 22 r/ 2 - 22400959860?? 2 v 26 + 26907093v 40 + 4639059v 44 + 126080307v 36 
+ 467481699v 32 - 25v 60 + 624309v 48 + 58263v 52 + v 64 + 2961v 56 - 2708022132v 34 ^ 
_ 740454048v 38 ?} 2 - 108490428v 42 r/ 2 - 164715 12v 46 r? 2 - 662580v 50 ^ 

- 9226708056v 30 t) 2 - 92880v 54 r/ 2 + 216v 62 r} 2 + 2916v 58 r} 2 + 2754990144 , 

(B.5) 

Q 6 = -12223143 - 20785248v 4 - 2196720v 24 - 6722352v 20 - 15363675v 16 

- 23024736v 12 - 23225940v 8 - 579312v 28 - 172186884t) 2 v 2 + 752520456^v 6 

- 108413964?? 2 v 10 + 166754376v 14 r) 2 + 257348664r) 2 v 18 + 140702832v 22 r? 2 

+ 37144008^v 26 - 828v 40 + 16v 44 - 14640v 36 - 114813v 32 - v 48 + 2349324v 34 ?7 2 
+ 425736v 38 r) 2 + 15876v 42 r? 2 + 648v 46 r) 2 + 12227760v 30 r} 2 , 

(B.6) 

Q 7 = 163046532875760v 32 - 108446607008910v 36 + 26922194450352v 40 

- 11404592789688v 44 + 3756782036502v 48 - 421486536163500v 28 

- 1062657813600v 52 - 45397293396v 60 + 290922851196v 56 - 1976417606251692v 20 
+ 2525883800478795v 16 - 2045958218223822v 12 - 68193367053291v 4 

+ 518822412814782v 8 - 1129718145924 + 703601320689252v 24 + 20v 84 - 2958074v 76 

- 1437073423v 68 + 13224991788v 64 + 175580978v 72 + 140423v 80 , 

(B.7) 

Q s = 8v 80 - 412v 7 6 + 1130201v 72 - 36420582v 68 + 661666116v 64 - 4050893934v 60 

+ 31716802512v 56 - 92949815898v 52 + 550493727984v 48 - 1678881463974v 44 

+ 4547507588514v 40 - 13314555136710v 36 + 9923931144288v 32 - 82854890580186v 28 

+ 1492061941800v 24 - 51956321216238v 20 - 427117486358664v 16 

+ 308048735088522v 12 - 243317414920563v 8 + 4198088418804v 4 + 1938652126956 . 

(B.8) 

In [15] an expression for f! 2 = c' 2 + \b' 2 was presented. It is straightforward to 
integrate it to obtain (imposing the proper boundary condition, i.e., /2(f) — > as 
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+ ): 



v 2 



. >/3(5i- 954) + , v 2 , 3475^3(^-3) 

f 2 = arctann^= H ; -r — arctan — — 

288 y/3 128(3 + v 4 ) y/3 



57o(108057o + 82 - 357q5i) 3 - r ' 



,.4 „,2 



In t—^ + tt^t, . K77 5 ^ -2125764^ 



24 3 + v 4 1152(3 + v 4 ) 4 {v 8 - 9) 4 

+ (14171765iC + 27209779277^! + 226521411840^ - 51018336^^2 + 408146688C 
+ 170061 125 3 )v 2 + (11809805iC + 272097792C + 158723712^^ _ 3968092877o<5 2 
+ 17006112(53 + 174686782464^)^6 + (11074698999 - 49601165i)w 4 + (-8468512335 

- 25981565i)t; 8 + (-3149285iC - 100147104^^ _ 2857026816C + 13226976^2 

- 816701522688^ - 1889568tf 3 )u 10 + (3564^ + 91909647)?/ 2 + (-165888^<Ji 

- 597756672^ + 37 32480C - 1512<SiC + 44064r} 5 2 - 285125 3 )?/ 4 + (866052^ 

+ 22334044221)*; 12 + (-74952864^^ - 4986365iC + 3355872768C + 17006112t7 <5 2 
+ 1203035037696^ - 69284165 3 )t; 14 + (6110429886 + 9972725i)t/ 6 + (419904^^2 
+ 83980807^1 - 1375605504C - 12597125 3 - 34992<5iC - 1340988618240^)t; 18 
+ (9797765 3 + 481310760960^ + 758165iC + 10077696C - 2659392?7 0( 5 2 
+ 1287705677^1 )v 22 + (874805i + 9442000710)v 20 + (62 2 8 007/ 2 - 2592t) 5 2 
+ 1728077^1 + 36<Ji(7 - 864<5 3 )?/ 8 + (-3125709 + Um x )v m + (-1440t7 O 5 2 + 8645 3 
+ 6324480^ + 485iC + 6624?y^i)i; 42 + (365i + 645597)?/ + (-12<Si + 105273)*; 
+ (480?7^i + 965 3 + 4(5iC - 967/0^2 + 5529607/ 2 )?/ 6 + (— 1108085i + 3107022786) 
+ (-1212246051847/ 2 - 419904t} (5 2 + 153964800C + 15552<5iC + 88 6 4 647/ 2 5i 
+ 3265925 3 )v 26 + (-25272^ + 208396314)*; 28 + (1710727) 5 2 - 917568?7^i 

- 42923520C + 7936745472t/ 2 - 45365iC - 466565 3 )t; 30 + 18648796131^ . 

(B.9) 



M 

)v 2A 



26 



We can further use a constraint at the second order (similar to (3.6)) to find 

10425^(^ + 2^ + 9)^ v 2 v 3 / ctV7cnA . 

4 = ^7-5 — txz arctan — + — — — — — — - 76527504<5 3 

2 8(v 8 -9) 2 72(w 8 -9) 5 (3 + v 4 ) 4 V 

+ 1836660096C + 1101996057600?^ + 994857552^ - 153055008r) 5 2 

+ 63772925iC + (-864^2 + 125iC + 576r$5i - 2016<5 3 )v 44 + (1190427845 3 

- 289103904r) 5 2 + 1851353376768r) 2 + 1938696768r) 2 ^ + 3265173504C 
+ 106288205iCo)w 4 + (-339560128512?? 2 - 3129840r} 2 5i + 674179200C 
+ 152285iC + 1905125 3 - 738720^ 2 )v 32 + (1912896t) 5 2 + 2034579326976?? 2 
+ 6065285 3 - 21384(5iC - 22581504?^i - 2058089472C )w 28 + (96435034542 
+ 17496005i)t; 22 + (194755671 - 17712(h)*; 38 - (24091992^ - 445385644605)t; 10 
+ (24383227336704?? 2 + 69838433280C - 56687045 3 + 2913084<?iC 

- 1795089607^2 + 1795089600r/^iK 2 + (-255091685i - 86716822293)w 6 

- lUS 3 v 48 + (18828186993 - 122821925!)*; 14 + (-9977140749312r/ 2 

- 78941952C - 71313696?} 2 5i + 16096320^2 - 2741045iC )v 24 + (5598720C 

- 86405 3 + 468^ C - 4752^^ - 28512r/ 5 2 - 618098688r} 2 )t; 40 + (314464150422 

- 16796165iy 8 + (-278645i + 41837067)v 3 4 - (57209092209 + 12754584<5i)t; 2 
+ (1788485i + 11908840482)t; 30 + (21413128771584r/ 2 - 3359232r} 2 5i - 54587525 3 
+ 18280940544C - 647352<?iC + 28133568r) 5 2 )w 20 + (7776<5 3 - 212544r/ 2 5! 

- 261792r/ <5 2 - 93312000Q, + 16858865664r/ 2 + 4644<5iC )t; 36 + (10380965i 
+ 32981460282)?/ 6 + (77944685iC - 33604077312C + 2406364848^<Ji 

- 8771888618496r} 2 - 300441312r/ 5 2 + 566870405 3 )w 8 + (-3365i - 357507)t; 46 
+ (-16(Ji - 4608)?/° + (-28973376^2 + 619148448^^ - 54268392960C 

- 154314725 3 - 36532030952448r} 2 - 524885iC )^ 16 + (-4094631 - 3384<5i)?/ 2 

(B.10) 

Using (3.2)-(3.6) we can obtain a decoupled equation for a 3 : 



5v 16 + I87; 12 + 216t; 8 + 126?/ + 243 
° " ° 3 + (^-9>(^ + 2^ + 3)(t; 4 -2^ + 3) a3 + J *° ' ( ' ; 
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u<13 ~ ~ (v 8 -9) 2 (v 8 + 2v± + 9) J2 ~~ 3(w 4 - 2v 2 + 3) 2 (v 4 + 2v 2 + 3) 2 (v 8 - 9) 6 (3 + w 4 ) 

, _ 16v 6 \/3^ 2 

X ° 3 ~ 9{v 8 - 9) 6 (3 + t; 4 ) 2 (t; 4 - 2v 2 + 3) 2 {v 4 + 2v 2 + 3 )2 arctanh ^ 

125100v 6 r?ov / 3^3 w 2 

arctan ■ 



(w 8 - 9) 4 (w 4 - 3) 2 (w 4 - 2v 2 + 3)(v A + 2v 2 + 3) ^ 

^4 , 3-t; 4 

27(v 8 - 9) 5 (v 4 - 3)(v 4 - 2v 2 + 3) 2 (v 4 + 2v 2 + 3) 2 H 3 + v 4 

5 



+ 
+ 



2fj 5iv 4 Tlr 



+ 



9(v 8 - 9) e (v 4 - 2v 2 + 3) 2 (v 4 + 2v 2 + 3) 2 

2v%TZ 6 

27(v 8 - 9) 5 (v 4 - 3)(v 4 - 2v 2 + 3)(v A + 2v 2 + 3) 

192v%5 3 n 7 
{v 8 - 9) 4 {v 4 - 3) 2 (v 4 - 2v 2 + 3){v 4 + 2v 2 + 3) 

4v%n 8 

(v 8 - 9) 7 (v 4 + 2v 2 + 3) 2 (v 4 - 2v 2 + 3) 2 (3 + w 4 ) 4 ' 



(B.12) 



with 



Tlx = 143489075i + 8264970432 + (-1435349865024 + 9087641 IS^v 4 + (462353675i 
+ 6473308508352)*; 8 + (-232062575i - 11722074916032)*; 12 + (11877204669696 
- 165927695iy 6 + (10431995i - 7723609200000)*; 20 + (164681l£i 
+ 3567784040064)*; 24 + (-852935i - 1321470640512)*; 28 + (—28431*1 
+ 440473837248)w 32 + (-132105017664 + 609935i)w 36 + (4293<5i + 31786872768)w 40 
+ (-5435190720 - 75875i)v 44 + (594864000 - 1179*i)v 48 + (2615i - 36391680)*; 52 
+ (933120 + 575i )v 56 + v m 5 1} 

(B.13) 
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TZ 2 = -16529940864r}o + 47829695 2 - 430467217/otfi + (138174665 2 + 532631427840t) 

- 149866362r} 5i)v 4 + (86802035 2 + 3526999604352^ - 101505231t) 5i)v 8 

+ (660871 1172096r}o - 16533725 2 + 77944687^1 K 2 + (-32476955 2 + 26394903r/ 5i 

- 1430009945856r} )^ 16 + (5786802^^ + 50962 10245632r) - 8004425 2 )t; 20 

- (1358127^1 + 277988205312^ - 2471315 2 )w 24 + (726514541568t} + 1807925 2 

- 542376?)o5i)t; 28 + (-30848 107392?) + 274595 2 - 13851r} 5i)w 32 + (-98825 2 
+ 62 8 92 1 1 3767}o - 12150r/ 5i)w 36 + (-1974761856??o - 4455<5 2 - 9477r) 5i)t; 40 
+ (-2525 2 + 324r) 5i + 1007023104r) )t; 44 + (1475 2 + 837^1 + 60466176r) )^ 48 
+ (1119744r}o + 126?^i + 265 2 )v 52 + (3^i + S 2 )v 5e , 

(B.14) 

U 3 = 7v 28 - 3v 24 - 93v 20 - 63v 16 + 1269v 12 - 4617w 8 + 6561w 4 + 2187 , (B.15) 

TZ 4 = -51656065200r) - 239148455 2 + (459165024^^ - 153055008?? 2 5 2 

- 165299408640?? > 6 + 71744535^1 + (596465856^ + 6845 2 - 2970^1 )v 44 
+ (-12397455648^ + 100442349r} 5i - 180689945 2 )t; 4 + (19486175 2 

- 3275759808720^0 + 39326634r} 5 1 )t; 8 + (-312232216320^ + 867311712^ 

- 289103904^5 2 y o + (-87^ - 102816r) - 34<5 2 )v 52 + (6022998r) 5i + 44877245 2 
+ 3920759121600?)o)^ 12 + (1539<5 2 - 3078?^i - 13050184752r} )w 40 + (28512^(5 2 

- 85536?7^i + 3079296077o> 42 + (1259712?)^ 2 + 1360488960r? 3 - 3 7791 36^ )v 22 
+ (67344203520^ - 18 706 72327^ + 62355744r? 2 5 2 y 4 + (368145^0^ - 128795 2 

- 282820705488r7o)^ 32 + (4115059207$ + 381024?? 2 5 2 - 1143072r7^i)w 38 - (115911*2 
+ 2547194665968^ + 2965572^1 )v 24 + (2592^^ - 933120?$ - 864?? 2 <5 2 )t; 50 - (55 2 
+ 15120r) )w 56 + (-641525 2 + 96228r} 5i + 143399735424r) )t; 28 + (50388480?? 3 

+ 46656r) 2 5 2 - 139968^5 1 )v 34 + (-7657^1 + 33<5 2 - 55887408t) )^ 48 + (440647? 3 5i 

- 15863040?$ - 14688?$<5 2 )?/ 6 + (26453952^! - 9523422720?$ - 8817984?$<5 2 )t; 26 
+ (92588832?$<5 2 - 277766496?$5i + 99995938560r) > 18 + (68283t/ 5i + 162<5 2 

+ 23803727904?) )?j 36 + (1924904811168t/ - 5570289^1 + 131225 2 )?j 20 
+ (8817984?$5i - 3174474240?$ - 2939328?$5 2 )?j 30 + (-112 1 931?^i 

- 9513261568080r}o + 11219315 2 )t; 16 , 

(B.16) 
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TZ 5 = -(7164612C + 4472292528r/ 2 )v 18 + 119062467v 16 + 223087122v 12 
+ 315675954v 8 + 74933181v 4 - 9198v 44 - 1583469v 32 - 82638v 40 - 430029v 36 

- 915v 48 + 13v 52 - 2032452v 28 + 34031907v 20 + 4452732v 24 + (115560C 

- 16708032?? 2 )v 34 + (-46530612C - 10489936752?? 2 > 10 + (17712C - 3045600r} 2 )v 38 
+ (52704?? 2 + 648C )v 46 + (12C + 1872r/ 2 )v 50 + (-118646208?? 2 + 470448C )v 30 

+ (IIOI6O77 2 + 2268C )v 42 + (-105 788 76487} 2 - 15903864C )v 14 + (-482772960?? 2 
+ 173016C )^ 26 - (25509168Co + 3520265184?? 2 > 2 - (2133532224?? 2 + 2799360C )^ 22 

- (127545840Q, + 21070572768r? 2 )v 6 + 23914845 , 

(B.17) 

U 6 = -885735 + (236196Q, + 62355744r} 2 )v 2 - 885735v 4 + (866052C 
+ 288159120?? 2 > 6 - 230291k; 8 + (-26453952^ - 157464C )v 10 - 177147v 12 
+ (66624768r? 2 - 134136C )v 14 + 275562v 16 + (-5412096r? 2 + 27216C )^ 18 
+ 92826v 20 + (9072C - 1594080r} 2 )v 22 + 7938v 24 - (4968C + 1684800r) 2 )v 26 
- 2214v 28 + (-31104?? 2 - 648C )v 30 - 963v 32 + (-50976?? 2 + 396C )v 34 - 531v 36 
+ (-1584?) 2 + 12C )v 38 - 3v 40 + v 44 , 

(B.18) 

n 7 = 7v 28 - 3v 24 - 93v 20 - 63v 16 + 1269v 12 - 4617v 8 + 6561v 4 + 2187 , (B.19) 
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n s = 1750802750915778v ib - 547909426568250v i2 + 551849751221706v 8 
- 3229924616793v 4 + 10220357329866v 44 + 252995054673834v 32 
+ 32427971352054v 40 + 94323062131416v 36 + 1924492298742v 48 + 491409014598v 52 
+ 392340630598470v 28 + 310374023297562v 20 + 903540325216554v 24 
+ 31603772430v 56 + 11854881954v 60 - 177976890v 64 + 55953441v 68 - 376v 76 
+ 185736v 72 + 16v 80 + (-77182514538916608^ - 361083701553408C )^ 14 
+ (6169681384704^0 + 202283833009397767? 2 )v 26 + (-39822513408r? 2 
+ 685117440Q> 62 + (117792358596864C + 19949787269471232t} 2 )v 10 
+ (211985675693568C + 64589051749911552r} 2 )v 18 - (23950080C + 403273728^) 
x v 66 - (4040529 154973568?? 2 + 29968782786432C )v 6 + (8159616^ + 746496C )v 70 
+ (24405013822464?? 2 - 7899545088C )v 50 + (107114016798720^ 
+ 297538935552C )^ 2 + (-281691 198720Q, + 348792330451968r? 2 )v 42 
+ (20993 178240C - 6791658865920r} 2 )v 54 + (-3560308992C + 487741077504t} 2 )v 58 
+ (-7095537792C - 244469231 1194112t} 2 )v 38 + (2463432321024C 
+ 3589000230500352r} 2 )v 34 + (-15313028838912C - 16874613168627456^)v 30 
+ (-1773530468928007} 2 + 190336884480C )^ 46 - (1384616551461 12C 
+ 66703936055747328r? 2 )v 22 + 697356880200 . 

(B.20) 

For the computation of C we actually need only the asymptotic solution of a' 3 (v) 
asx = 3 1/4 - v -> 0+. Using (B.11)-(B.20) we find: 




(B.21) 



where A3 must be fixed so to satisfy the boundary condition (3.1). 
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C Source terms for the framework of Kinoshita et al 



Order - 0: 

S° A = 180y w (16 - 23y 4 ) , 
S° B = -432(h/ 10 , 

S° c = 0. (C.l) 

Order - 1: 

S\ = -(sSw^ + y' + y + iyy 1 

% 6 (1L4° (y 3 + 2y 2 + 3y + 4) - 3y 4 (74945y n + 191854?/ 10 + 308763y 9 

+425672/ + 370056y 7 + 209256y 6 + 48456y 5 - 112344?/ 4 - 153504y 3 

-109728y 2 - 65952y - 22176)) , 
! 8640|/ 10 (ll + 26y) 

^ = -(SS^y-^^ + ^ + y + l) 3 )- 1 

(2y 3 (3 (84y 10 + lQ8y 9 + 252y 8 + 2541y 7 - 754y 6 - 4049y 5 - 7344y 4 
-15624y 3 - 12528y 2 - 9432y - 6336) y 7 + 11 A$ (3y 2 + 2y + l))) . (C.2) 
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